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Abstract 

o ■ 

' The trace over the degrees of freedom located in a subset of the space transforms the vacuum 

state into a mixed density matrix with non zero entropy. This is usually called entanglement entropy, 
and it is known to be divergent in quantum field theory (QFT). However, it is possible to define 
a finite quantity F(A, B) for two given different subsets A and B which measures the degree of 
entanglement between their respective degrees of freedom. We show that the function F(A, B) is 
I . severely constrained by the Poincare symmetry and the mathematical properties of the entropy. In 

particular, for one component sets in two dimensional conformal field theories its general form is 
completely determined. Moreover, it allows to prove an alternative entropic version of the c-theorem 
for 1 + 1 dimensional QFT. We propose this well defined quantity as the meaningfull entanglement 
entropy and comment on possible applications in QFT and the black hole evaporation problem. 

7— I ! 

1 Introduction 

-xt- : 

The density matrix |^')(^ f | corresponding to the vacuum state \^>) can be traced over the degrees of 
■ freedom located in a subset of the space giving place to a mixed density matrix with non zero entropy. 
This is usually called entanglement or geometric entropy. It was argued that this quantity is deeply 
related to the entropy of black holes [TJ. In the black hole spacetime the causal structure provides a 
• natural partition of degrees of freedom, which corresponds to the asymptotic data at infinity and at 
J> , the horizon respectively. The quantum state relevant for asymptotic observers is obtained by tracing 
k>( | the complete density matrix over the invisible degrees of freedom on the horizon. But contains 
correlations between the field modes inside and outside the black hole, and in consequence, the partial 



trace leads to a mixed density matrix with non zero entropy. Indeed, several different calculations in 
specific examples have shown that the geometric entropy is proportional to the area in Minkowski space 
[H El E| • More recently the entanglement entropy has been also a subject of much interest in the 
study of condensed matter systems, specially in relation with the density matrix renormalization group 
method [6] . 

More generally, the entanglement entropy of the vacuum state can be thought as a special case 
of a more familiar quantity, the entropy S(V) contained inside a volume corresponding to a given 
region V of the space. For gases and other systems this is usually calculated by defining the system 
in question in the bounded region, what requires the specification of boundary conditions. However, 
these normally kill the correlations that the global state p could present between degrees of freedom 
in V and its complement —V. To compute S(V) without artificially modifying the global state p 
we need a partition of the total Hilbert space H as a tensor product of Hilbert spaces inside and 
outside V, Hy H_y. Then, the local density matrix is defined by py = Trn_ v p and its entropy is 
S(V) = — Trpy log py. The local Hilbert spaces Hy, density matrices py and entropies S(V) must 
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Figure 1: (a) The spatial surfaces A and A' have the same causal domain of dependence given by the 
diamond shaped set in this picture. The same entropy should correspond to A and A'. Thus, we do not 
make distinctions between a diamond set and any of its Cauchy surfaces, (b) Two spatial surfaces A 
and B included in the same global Cauchy surface. Writting area(X) for the d — 1 dimensional volume 
of the boundary of a d dimensional spatial set X we have area(^4)+area(i?) = area(^4ri-B) + area(^4U-B). 

satisfy certain compatibility conditions that simply express that they arise from a global state p and 
Hilbert space H by partial tracing. Given two non intersecting sets A and B we have that 

Haub = Ha®H b , PA = Tr HB p AuB ■ (1) 

These conditions allow us to use a very important mathematical property of the quantum entropy 
which is called the strong subadditive inequality (SSA) and which is usefull to relate the entropies 
corresponding to the different subsystems involved in a tensor product [8]. In the present case it 
acquires the following geometrical form [U [9j [10] 

S(A) + S{B) > S{A n B) + S{A U B) , (2) 

for any two sets A and B. In addition, the purity of the vacuum state gives the duality relation 

S(A) = S(-A). (3) 

In Minkowski space the independent degrees of freedom should be assigned to subsets of general 
Cauchy surfaces and an entanglement entropy must correspond to any of these subsets. Moreover, 
the causal structure and unitarity imply that this entropy must be the same for different spatial sets 
having the same causal domain of dependence (see fig. 1(a)). It has been shown that the combination 
of the SSA property, the positivity of the entropy, and the Poincare invariance of the vacuum state in 
Minkowski space strongly constrain the function S in any dimensions |10j . The result is that if S is 
finite for at least one arbitrary set then it has to finite and proportional to the boundary area for any 
other set, plus a constant term: 

S(A) = s avea(A) + a . (4) 

Here area(^4) means the area of the boundary of the spatial surface A and s and a are positive constants. 
Thus, the area law for the entanglement entropy has its origin in very basic and fundamental properties 
of relativity and quantum mechanics. In 1 + 1 dimensions the equation corresponding to (j4]) reads 

S(A) = s niA + a , (5) 

where uia is the number of connected components (or the number of boundary points) of the one 
dimensional set A. The SSA relation was originally introduced to show the existence of the mean 



2 



entropy (the limit of entropy divided by volume for large sets) in translational invariant systems [U [9] . 
We stress that the presence of additional boots symmetries makes the SSA relation a much more 
powerfull tool. For example, the equation ([5]) tells that the entropy of a single component set in 1 + 1 
dimensions is a constant (if it is finite), while the SSA relation combined with translational symmetry 
alone only implies that this is a concave function of the length. For different types of sets or more 
dimensions the contrast is much greater. 

The equations and ([5]) apply for finite S. However, all specific calculations, numerical and 
analytical, have shown that the entanglement entropy is ultraviolet divergent, even for free fields. This 
divergences can not be softened by adding different matter fields since the entropy is a positive quantity. 
This indicates that the entanglement entropy may not have a finite covariant meaning in quantum field 
theory (QFT). 

The calculations of S present in the literature have shown that the divergent terms are proportional 
to the area, which is also suggested by the theorem leading to eq. @) [TO]. Motivated by this fact, 
we introduce a new quantity constructed out of the function S which is however free from divergences. 
Let A and B be two spatial sets in Minkowski space which belong to the same Cauchy surface QI (see 
fig. 1(b)). The symmetric function 

F(A, B) = S(A) + S(B) - S(A n B) - S(A U B) (6) 
is finite, since the divergences cancel due to the relation 

area(A) + area(-B) = area(^ U B) + area(yl n B) . (7) 

It is also positive because of strong subadditivity and satisfies the duality 

F(A,B) = F(-A,-B). (8) 

Besides, F(A, B) has a very suggestive property, not shared by S. When A and B are non inter- 
secting sets F(A,B) is monotonically increasing, 

F(A,B) < F(A,C) , BCC , CPiA = 0. (9) 

From all this, we are lead to propose F as a well defined measure of the degree of entanglement in 
QFT. In fact, for two non intersecting sets A and B it can be written as 

F(A,B) = S{A) + S(B) - S{AuB) 

= Tr(p AuB logpAuB) ~ Tr(p AuB log(p A <8> Pb)) = S(p A \jB \ PA ® Pb), (10) 

where S(pi \ p?) = Tr (pi (log (pi) — log(p2)) is known as the relative entropy for two states p\ and p2 
acting on the same Hilbert space. This later can be thought as a measure of the statistical distance 
between a pair of states. In the context of quantum information theory, the particular relative entropy in 
(I10p is called the mutual information between the states p A and ps in the composite quantum system, 
and it is used as a measure of the amount of information they have in common |12j . The function 
F(A, B) for intersecting sets can be defined in terms of the relative entropy for non intersecting ones 
with the use of additional auxiliary sets. 

The connection between F and the relative entropy is also satisfying from the mathematical point 
of view. The origin for the divergence of the standard entanglement entropy can be traced to the 
impossibility of expressing the Hilbert space H as a tensor product of Hy <X> H—y in the relativistic 
case. This partition is unambiguous if the system is defined on a lattice, the local Hilbert spaces 
being generated by the local degrees of freedom (for example the spin- like operators in a given region). 

1 Constructions with two sets where also used in a related context to calculate quantum fluctuations of operators 
attached to the local regions [11] . 
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However, in a relativistic QFT the axiomatic investigations have shown that although there is a well 
defined notion of local algebras of operators in a volume V (type III Von Neumann algebras in general) 
there is no Lorentz invariant partition of the total Hilbert space into a tensor product Hy ® H_y (see 
|13} I14j). One heuristic reason for this can be seen considering the theory of a free scalar field 4>. In 
the classical theory the phase space attached to a bounded region V is given by initial data <j)(x) and 
4>(x) vanishing outside V. In the quantum case however, only positive energy solutions of the wave 
equation appear in the annihilation part of the operator 4>{x). In momentum space the time derivative 
is then e p = —i\/p 2 + m 2 . Now, in real space, the operator given by e p is antilocal, what means that it 
transforms any function vanishing outside an open set V in a function which does not vanish identically 
in any open subset outside V [2j. This implies that one can not have (0\4>(x)\ip) and {0\4>(x)\ifj) 
vanishing simultaneously in the same region of space, where is any one particle state. In general 
the moral is that in QFT there is no covariant meaning for the localization of states, and what can be 
unambiguously localized are the field operator algebras. In this sense it is remarkable that even if a 
notion of entropy for a state in a given algebra may not exist (type III algebras do not contain a trace), 
the relative entropy of two states for the same Von Neumann algebra has a well defined mathematical 
meaning [15J. 

2 Entanglement entropy in 1 + 1 dimensions 

The equation © implies the following relation between the finite F functions 

F{A,B) + F(A,C)+F{AuB,AuC) + F{AnC,AnB) 

= F(B,C) + F(A,BuC) + F(A,BnC), (11) 

for any A, B and C. In a given regularization scheme, as S goes to infinity these relations between finite 
quantities are maintained. The most general form of a function that solves (llip in 1 + 1 dimensions 
can be found by a long but straightforward calculation that we are not presenting here. It follows that 

F(A, B) = G(A) + G(B) - G(A n B) - G(A U B) . (12) 

Evidently, this equation has the same form as ([6]) but now the one set function G is perfectly finite. 
Since the meaningfull quantity is F, the function G(A) is defined up to an arbitrary term proportional 
to the number tua of connected components in A. Thus, for the one component sets there is a "gauge" 
symmetry associated to an additive constant. The positivity of F implies strong subadditivity for G, 

G(A) + G(B)>G(ADB) + G(AUB), (13) 

and from (JSJ) we also have G(A) = G{-A). 

The crucial difference between G and an entropy is that it can attain negative values. This must be 
so if G has some dependence on the size of the sets, because otherwise, as we have already mentioned, 
eq. ([5]) would apply for G. Moreover, it must attain negative values with arbitrarily big modulus. If 
not, it could be converted into a positive function by adding a suitable term proportional to txia- 

2.1 Entanglement entropy in Conformal Field Theory 

For the special case of a conformal field theory (CFT) and when the sets A and B are one component 
and intersecting, the eq. (|12j) allows us to obtain the general form of F(A, B). We see from figure 2(a) 
that only four points determine the position of the diamonds corresponding to A and B. As F(A,B) 
must be invariant under global conformal transformations, general results ensure that it must be a 
function of the cross ratios 

F(A,B) = F(r h ,r h ). (14) 
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Figure 2: (a) Two intersecting one component sets A and B whose straight Cauchy surfaces have sizes 
and tb respectively. The diamonds formed by the intersection and union (followed by causal completion) 
of A and B have sizes tav\B and taub- This configuration of sets is uniquely determined by the position 
of the points x\,X2,x% and X4 forming the spatial corners of the diamonds, (b) Configuration of two 
intersecting sets of size ^Jr\r<i from which eq. (|19p can be obtained by means of the SSA inequality. 



with 

^23 ^14 v 23 v 14 
Vu = , Vv = , (15) 

U13 U 24 V 13 V 2 4 

where we have written u%j = u% — Uj and Vij = v% — Vj and the (iii,Vi) are the null coordinates of the 
point X{. On the other hand, for one component sets, G can be written as a function of the length 
of the diamond base = y 7 (xi — Xj) 2 = ^/WjUij- Using this and equations (fT2|) and (fl"4l) we get the 
general form for G and F 

G{A) = k\og{r A ) + f3. (16) 
F(A,B) = felogf TAVB ) , (17) 

\TAnBrAUBj 

where k and (3 are constants, with k > 0. 

The conformal symmetry provides a simple argument in order to explain why F is finite while S is 
not. Both quantities are supposed to be invariant under global conformal transformations which leave 
the vacuum invariant, but while F is a function of four points only two points determine the argument 
of S for one component sets, and there is no possibility of forming a cross ratio with two points. Thus 
S must be either infinite or a constant. 

Several authors have calculated the entanglement entropy S in 1 + 1 conformal field theories with 
different regularization prescriptions, both numerically and analytically [21 HI El HE]. The result is 
S = (1/6) (c + c)log(r^/e), where e is an ultraviolet cutoff, and c and c are the holomorphic and 
antiholomorphic central charges. This formula can not be used for I a < e that gives S < 0. According 
to equation ([6]), our results (|16p and (|17p . obtained from very general considerations, are perfectly 
consistent with this expression for S, identifying k with (c + c)/6. Thus, G can be interpreted as a 
regularized entropy, G(A) = fclog(r^), with fi a renormalization scale. The fact that G is negative in 
some range causes no problem since its physical meaning is only up to a constant term. 

From eq. (|17p we see that the divergences appear in F reduced to certain special configurations 
of the diamonds. This is when r^nB goes to zero, which is not surprising, since in that limit the two 
diamonds separate each other and F turns into a different function, the one for non intersecting one 
component diamonds we discuss below. 
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2.2 An entropic c-theorem 

Now we look for the general form of G(A) = G(ta) for one component sets in 1 + 1 dimensional 
theories. Consider two one component sets A and B of the same size r% , intersecting in a segment of 
size r\ — r 2 with r 2 < r±. It follows from (|13|) applied to A and —B that 



G(n)>G(r 2 ) , n>r 2 . (18) 

Thus, G(r) is increasing and G' > 0. The SSA relation for intersecting diamonds disposed as shown in 
the figure 2(b) reads 

G(VnT-2)>l(G( ri ) + G(r 2 )). (19) 

It can be shown that all other inequalities coming from SSA for two intersecting one component dia- 
monds are less constraining than this relation. The meaning of (|19p becomes more clear defining the 
function 

G(a) = G(e a ) , (20) 



such that the domain of G is the whole real line. Equation (|19|) becomes 

G{\{a + b))> l -(G(a) + G(b)). (21) 

This relation is called weak concavity, and it is known to be a property of the entropy for translational 
invariant systems in the original length variable. Here we have it in the logarithmic variables. Using 
repeatedly this equation we can prove concavity for G, that is 

G(\a + (1 - X)b) > \G(a) + (1 - \)G(b) , (22) 

for any A £ [0, 1]. The concavity of G can be more simply written G" < 0. Introducing a new function 
C(r) = r G' (r) we summarize all the constrains on the original function G as 

C(r) = rG'(r) > 0, (23) 
C'{r) = G'{r) + rG"{r) <0. (24) 

According to the equations ()23[) and (|24|) the dimensionless function C must be positive and de- 
creasing with r. Since it depends on G' it is also a physically meaningfull quantity independent of the 
ultraviolet cutoff. Thus, the renormalization group equation for C reads 

i 

where r is the substraction point, the g$ are the different dimensionless coupling constants and Pi(g) 
the corresponding beta functions. This equation expresses that the total derivative of C with respect to 
t is zero, what follows from the fact that C is an observable quantity in a renormalizable theory. Note 
that in consequence, the anomalous dimension term is not present in (|25p . By dimensional analysis C 
has to be a funtion of (rr) and the dimensionless coupling constants, then we have 

r *-^) c =°- (26 > 

As C is decreasing with r we see from eq. (|25p and (|26p that C(ro,g) for fixed ro is a function of 
the coupling constants which decreases along the renormalization group trajectory. Moreover, at the 
fixed points, corresponding to conformal field theories, eq. (|16p implies that C(ro,g) is a constant 
k = (c + c)/6, stationary under the renormalization group flow. This is in fact the statement of the 
Zamolodchikov theorem 1171. 
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Remarkably, this theorem follows here from the properties of the entropy and the Poincare symmetry 
without resorting to quantum fields. Also, the unitarity and causality of the theory are implicitly used 
in the fact that the entropy of a given set does not depend on the particular choice of the Cauchy 
surface, what is crucial for the arguments that lead to eq. (|19p (see figure 2(b)). Moreover, the SSA 
property requires a self-adjoint and positive density matrix which is also related to unitarity (however, 
it also works for classical probability distributions) [8]. 

The C charge in eq. f|23[) has an interesting physical meaning: it measures the variation of the 
geometric entropy with the size of the set whose degrees of freedom are traced over. For example, for 
massive fields, the entropy saturates at a given radius (see for example [5]), what leads to the vanishing 
of C for large r. This is in accord with the renormalization group flow interpretation of the running of 
C given by eq. (|26p since at the infrared no available degrees of freedom remain in the massive case. 

In relation with the massive C function we make the following remark that gives additional support 
to the connection of the entropy F with the algebraic approach to QFT. Numerically we find that the 
C function ()23f) for a free scalar field of mass M decreases exponentially, C(r) ~ e —V3Mr^ f or } ar g e 
r, what also happens for the Zamolodchikov c function. Adding independent scalar fields we have a 
contribution to C(r) due to the sum of the exponential tails which is of the form J dMa(M)e~^ Mr , 
where a(M) is the density of fields per unit mass. We can choose an exponentially increasing density 



cr(M) ~ e T H in such a way that the theory has a maximal temperature Th (Hagedorn temperature). 
In this case the divergence of G does not take place at r = but at a finite distance r ~ TZ . This 
coincides with the results of [18] which tell that a consistent tensor product structure can not be given 
to the states on the algebras of two non intersecting sets when their separation distance is less than 
some constant times T7r . 

Some relations between the geometric entropy and the renormalization group flow were already 
suggested in the literature based on numerical results [5]. Moreover, interesting entropic c-theorems 
were introduced in [19]. These also rely on the properties of the entropy. However, they do not involve 
the entanglement entropy between the degrees of freedom in different regions of the space but the 
relative entropy between the vacuum states corresponding to different values of the theory couplings. 

The renormalization group flow of the Zamolodchikov c function is given in terms of the energy 
momentum tensor correlators. It is likely that the energy momentum tensor is also involved in the 
running of the entropic function C proposed in this paper. This would be very interesting, specially 
taking into account that one of the originals aims of the investigations about the geometric entropy 
was the underestanding of the relation between entropy and gravity. 

2.3 The two component set function 

The equations (|23p and (|24p give the most general form of the function F(A, B) in 1 + 1 dimensions 
for intersecting one component sets A and B. Of course, these leave a lot of freedom, but this must be 
so, since F contains information on the mass spectrum of the theory. The analysis can be continued 
considering the multicomponent set functions, that is, the case in which the sets A and B are multi- 
component or non intersecting. Then, the function F and the relevant inequalities involve the function 
G evaluated on multicomponent sets. The general analysis is outside of the scope of this letter. Here 
we show only some constrains for the simple case of two non intersecting one component sets A and B 
lying in a single spatial line in a CFT. These sets are as before determined by the position of their four 
edge points, what implies due to global conformal invariance that F is a function of the cross ratio 



M 



F(A, B) = F{rj) = k \og(r A ) + k log(r B ) - G(A U B) 



(27) 



with 



r A r B 



(28) 



V = 



(r A + r c )(r c + r B ) 
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where and rs are the sizes of A and B and tq is their separating distance. We can rewrite eq. ([2 
as 

G(A U B) = k log (— ; — - — -- — - — _ 29) 

V ( r A + r c)\ r C + r B) / 
with U(rf) = F{rf) + fclog(l - rj). The duality © leads to 

U(rj) = U(l - rj) . (30) 

Besides, from the SSA relation between AuB and another set C with one or two components it follows 
that 

U(rf) > fclogfa), (31) 

U'(V) < (32) 
V 

TT"( \ ^ fc /ooN 

17 w - -^t-^ki^f (33) 

These inequalities highly constrain the function U, and in particular show that it has to be bounded in 
its domain r\ E [0, 1]. Thus, the only singular points in G(A U B) come from the first term in the right 
hand side of (J2SJ) - In particular in the limit B — ► —A, tq = e — > we get F — ► 2/c log(r^/e). From eq. 
([6]) this gives 

5(r A )^A;log(/ A /e), (34) 

consistent with previous calculations [16]. 

We have checked numerically the formula (|29p for a Majorana and a Dirac fermion (c = 1/2 and 
c = 1 respectively) using the methods of [5j [20]. The result is f7 = and k = 1/6 and fc = 1/3 
respectively. On the other hand, we have also considered the theory of a free massless boson (c = 1). 
The existence of a zero mode provides an additional infrared divergence to the entanglement entropy 
which is also present in F (this happens for all theories with non compact global symmetry groups). To 
regularize it we have introduced a small mass M and attached the theory to a cylinder of circumference 
L such that ML <C 1. In this case the numerical result is 

U( Vcyl ) = i log(7 7cyl (l - r/ cy i)) - l - log(ML), (35) 



where 

(36) 



sin(^)sin(^) 
sin( 7r(rA + rc) )sin( 7r(rfl + ro 



is the expression of r\ in the cylinder. With r& ,^b , r c C£ one recovers r/ cy i — > r\. The formula (|35p 
satisfies the inequalities (I31II33P with = 1/3, except when rj < (ML)' 9 / 2 ' or (1 — rj) < (ML)' 9 / 2 ' where 
the mass enters into the game and equation (|35j) and the conformal symmetry do not hold any more. 
This is why (|35p does not seem to lead to a bounded U. The point we want to make here is that the 
function U has the ability to distinguish between different models with the same central charge (c > 1). 



3 Final remarks 



One could wonder if the whole set of entanglement entropies are in one to one correspondence with the 
different QFT. In other words, is it possible to uniquely define a QFT by the set of its entanglement 
entropies?. This is likely, since the local algebras of all QFT are the same mathematical objects, 
independent of the theory and the space-time dimension |13j . Thus, the way they are entangled to 
each other must play a central role. If one accepts this, a second question naturally appears: is it the 
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Figure 3: The Rindler wedge A and an exterior diamond set B. The curves of constant v are lines 
through the origin. 

only requisite of a consistent set of entropy functions to obey the SSA relations and the space-time 
symmetries?. 

Further interest to these questions is given by the fact that the entanglement entropy F has a nice 
property suited to study dimensional reduction. Take for example the sets of the form A x [x] in d + 1 
spacial dimensions, where [x] is an interval of length x and A is any d dimensional set. Then the 
function 

- . , , , F(A x \x\,B x \x\) , s 

F(A, B) = lim — i ^ ^ (37) 

x— »oo x 

is well defined and must satisfy all requirements of an F function in d dimensions [10] . Generalizations of 
the results of this paper to higher dimensions using this construction are currently under investigation. 

Finally, we suggest an interpretation of F in the black hole evaporation problem. Take A as the 
Rindler wedge and let B be a diamond outside A as shown in fig. 3. Then F(A, B) can be interpreted, 
at least in certain limit, as the amount of entropy from the Unruh thermal radiation that an observer 
can detect in an experiment running for a finite period of time corresponding to the set B (see related 
ideas in [HUT]). This period measured in the boost parameter v = arcth(i/x) is 

A„ = lo g (^£). (38) 

Now, the equations (|271) and (|29l) for ta — ► oo give 

F(A, B) = kAv + U ( — ^ — J . (39) 
\r B +r c J 

Thus, as U is bounded, for big enough Av this gives F(A, B) ~ kAv. We also have Av = aAr, where 
r is the proper time meassured by an accelerated observer with constant acceleration a. This formula 
tells that the entropy is proportional to the diamond time At times a constant flux of entropy ka due 
to a constant flux of Unruh radiation. Indeed, the Unruh temperature is T = a/(2n) and in the high 
temperature limit (adequate to the limit taken above) the flux of entropy per unit time in a CFT has 
a universal expression given by As/ At = tt^^T [22] . Therefore, the identification of F with the total 
entropy flux during the interval At lead us again to the equation k = In the black hole case Av 
should be interpreted as At/(4Mb#), where r is the time as measured by the asymptotic observers 
and Mbh is the black hole mass. Thus, F would be proportional to the asymptotic time, corresponding 
to a constant flux of Hawking radiation. For a real black hole the back-reaction must be included. In 
four dimensions the black hole lifetime is At ~ Mbh' giving place to a maximal observable entropy 
i^^-B) ~ M BH , which is proportional to the horizon area. 



9 



4 Acknowledgments 



H.C. thanks the invitation to visit the ICTP where this work has been done. 

Note added - After the submission of this paper to the arXiv data-base, the reference [23] appeared, 
showing a calculation of the multi-set CFT entropy functions. Their general results are in agreement 
with the particular cases considered in Sections 2. A and 2.C. 

Erratum: 

1- Equation ([7]) eliminates the divergences in F{A,B) only in two dimensions. In more dimensions 

additional subleading terms are present and F(A,B) does not converge for generic A and B. 

2 - In Section 2.B the long distance behaviour C(r) ~ e~^ Mr should read C(r) ~ e~ 2Mr . 

3.- In Section 3 the large Av limit is the low (rather than high) temperature limit. There is a factor 

two in As /At which should read ir^^T. This means that F represents twice the total entropy flux 

during the time At. 
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